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Possible ordered states in the 2D extended Hubbard model with on-site ([/ > 0) and nearest- 
neighbor {V) interaction are examined near half filling, with emphasis on the effect of finite V . 
First, the phase diagram at absolute zero is determined in the mean field approximation. For 
< 0, a state where d3,2_y2-wave superconductivity (dSC), commensurate spin-density-wave 
(SDW) and 7r-triplet pair coexist is seen to be stabilized. Here, the importance of 7r-triplet pair 
on the coexistence of dSC and SDW is indicated. This coexistent state is hampered by the phase 
separation (PS), which is generally expected to occur in the presence of finite-range attractive 
interaction, but survives. For V > 0, a state where commensurate charge-density- wave (CDW), 
SDW and ferromagnetism (FM) coexist is seen to be stabilized. Here, the importance of FM 
on the coexistence of CDW and SDW is indicated. Next, in order to examine the effects of 
fluctuation on each mean fleld ordered state, the renormalization group method for the special 
case that the Fermi level lies just on the saddle points, (7r,0) and (0,7r), is applied. The crucial 
difference from the mean field result is that superconductivity can arise even for U > Q and 
V > 0, where the superconducting gap symmetry is ci^2_y2-wave for U > iV and s-wave for 
U < 4:V. Finally, the possibilities that the mean field coexistent states survive in the presence 
of fluctuation are discussed. 



KEYWORDS: 2D extended Hubbard model, coexistent state, d-wave superconductivity, s-wave superconductiv- 
ity, spin-density wave, charge-density wave, ferromagnetism, vr-triplet pair, rj-singlet pair, phase 
separation 



§1. Introduction 

In connection with the studies of the copper oxide 
high-Tc superconductors with Cu02 planes, the elec- 
tronic states in two-dimensional systems has been in- 
tensively studied. Especially, the possibility of various 
ordered states has been discussed. One characteristic 
feature is the proximity of superconductivity and an- 
tiferromagnetispp. In our previous work (hereafter re- 
ferred to as I), El we have shown in the mean field ap- 
proximation that the coexistent state with d-wave su- 
perconductivity (dSC), commensurate spin-density- wave 
(SDW) and 7r-triplet pair can be stabilized near half fill- 
ing by repulsive backward scattering ('Umklapp' and 'ex- 
change') processes between electrons around the saddle 
points (tTjO) and (0,7r). As we shall show later, this model 
with such a particular type of interaction have similar 
features to those in a square lattice model with on-site 
repulsion U > and nearest-neighbor attraction V < 0, 
i.e., an extended Hubbard model.a Therefore, it is inter- 
esting to examine in more detail the possibility of the 
above coexistent state by use of this model. At the same 
time, the extended Hubbard model with both on-site and 
nearest-neighbor repulsion {U > and y > 0), is also 
of physical interest. In the 2D extended Hubbard model 
ioi U > 0, it has been shown based on the mean field 
approximation that extended-s-, p- and d-wave super- 
conductivity can arise depending on the electron density 



n for V < 0,BI3'dI3-' and commensurate charge- and spin- 
density wave (CDW and SDW) can appear at half filling 
n = 1 for V > O.clcP However, the property of the ground 
state for finite carrier doping and the relationship among 
various order parameters, especially between dSC and 
SDW, have not been understood yet, even in the mean 
field approximation. From these points of view, we will 
study possible ordered states, especially possible coex- 
istence of different orders, in the 2D extended Hubbard 
model on a square lattice near half filling for U > and 
V 0, with emphasis on electrons around the saddle 
points (7r,0) and (0,7r). 

In the extended Hubbard model is introduced and 
its relationship to our previous model used in I is referred 
to. In §|[ the phase diagram at absolute zero, T = 0, is 
determined in the mean field approximation. In §^ the 
effects of fiuctuation on the mean field ordered states are 
examined based on the renormalization method applica- 
ble only for the special case that the saddle points (7r,0) 
and (0,7r) lie just on the Fermi surface. 

§2. Extended Hubbard Hamiltonian 

The extended Hubbard Hamiltonian, H = Hq + Hjj + 
Hy, is written as follows: 

Ho ^J2^P^P'^^P'^' (2-1^) 



Hu = U^n^'^n,i^^^nq^n^ql, (2.1b) 
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-qi 



(2.1c) 



where a is the spin index taking a value of +1 (—1) 
for t (i) spin, and the opposite spin to a is denoted by 
W = —a. N is the total number of lattice sites, — 
ep — ^ is the one-particle energy dispersion relative to 
the chemical potential /i, including nearest-neighbor- (t) 
and next-nearest-neighbor- {t') hopping integrals, 

€p — —2t{cospx + cospy) — 4t' cos cospy, (2. Id) 



Efccr cLcfc-i-gT, P = ±S:,±y is the unit 



lattice vector and 



Vq = V{cosqx 



(2.1e) 



The energy dispersion Cp has two independent saddle 
points, (tTjO) and (0, tt). In this paper, weGxt'/t = — 1/5 
with t > 0, in which case the Fermi surface in the ab- 
sence of interaction approaches (tt, 0) and (0, tt) as the 
hole doping rate, (5 = 1 — n, is increased from half filling, 
6 = 0. 

Here, we examine the relationship between the 'g- 
ology' model used in I and the present extended Hubbard 
modelJa' In I, we have treated the backward scattering 
with large momentum transfer between electrons around 
(7r,0) and (0, tt), i.e., 'Umklapp' {g3±) and 'exchange' 
igi±) processes, and considered three types of the scat- 
tering processes, i.e., (II Cooper-pair, (2) density- wave 
and (3) 7r-pair channels. EP (Here we denote the Hamilto- 
nian for these channels as Hi, H2 and H^.) Especially 
for the repulsive case, gzi_ > and gi± > 0, we have 
shown that the coexistent state with dSC, SDW and vr- 
triplet pair can be stabilized near half filling at low tem- 
perature. However, the above eflFective interaction is too 
simplified in that (a) forward scattering processes are not 
included and (b) the fc-dependence of interaction is ig- 
nored. Therefore, by transforming Hi {i = 1,2,3), into 
real-space representation, we will obtain a well-defined 
model on a square lattice. If we keep on-site and nearest- 
neighbor density-density Coulomb interaction, we obtain 



Hi 



93±_ 

2N 



9± \ 



<ij> <y 



H. 



9i±_ 
2N 



EE^ 



(2.2b) 



(2.2c) 



<ij> a 

where (7j_ = 33 j_ -I- g\^, < ij > stands for a bond con- 
necting site i and its nearest-neighbor site j and a = 
16/7r* ~ 0.164. Each Hi is seen to describe on-site re- 
pulsion and nearest-neighbor attraction for g3±,gi± > 0. 
Therefore, also in the extended Hubbard model for U > 
and V < 0, the coexistent state with dSC, SDW and vr- 
triplet pair is expected to be stabilized. 

§3. Mean Field Analysis 

First, we determine the phase diagram at absolute 
zero, T = 0, near half filling in the mean field approxima- 



tion. We fix the electron density to n = 0.9. With this 
choice of parameters, the Fermi surface of noninteracting 
electrons is of the YBCO- or BSCCO-type and lies close 
to the saddle points (7r,0) and (0,7r), as shown in Fig. 1. 
t'/t=-1/5, n=0.9 



Fig. 1. The Fermi surface in the absence of interaction for t' /t - 
-1/5 and n = 0.9. 



3.1 Nearest- Neighbor Attraction < 

We start with the case V < 0. As we saw in §||, the or- 
dered states with dSC, SDW(= antiferromagnetism, AF) 
and TT-triplet pair is expected. We consider the following 
order parameters, 



CiaCi^p_(j ^ — 



(3.1a) 

amcos{QRi), (3.1b) 



"■Sp + 9p cos{QR. 
11 
2 



where Q = (tTjTt), s_p = Sp and g_p = qp. 



(Qp) 



stands for a spin-singlet (triplet) pair of two electrons 
with a total momentum (Q) and total spin S = 
(S* = 1 and Sz = 0), i.e., Cooper-pair (7r-triplet pair), 
and TO for the local staggered spin moment. We take 
Sx = —Sy = sq = real and = —q-g = qo = real, 
i.e., consider dj.2_y2wave pairing, which is favored near 
half fiUing.B While sq describes dSC, jOs, describes an 
electron-pair of Px — Py-wave symmetry.LptP This can be 
easily seen by writing the operator of 7r-triplet pair in 
fc-space. 



Ott — fr: ^ ' ^ ^ WpC—p+QaCpc , 



V2 



1 



pa 



= -7^Y.'^%+p%-paC%+p,. 
^ pa 



(3.2a) 



(3.2b) 



where Wp cx cospa: — cospy is the (ij.2_j,2-wave factor 
and < >oc go- The orbital function as a function 
of the relative momentum n of two electrons, wq , , is 
proportional to sinp^; — sinp^. 

The order parameter in the mean filed Hamiltonian 
are given as follows in terms of sq, qo and to, 

Adsc = -2\V\.so, A, = -2\V\qo, (3.3a) 
AsDW = Um, (3.3b) 

where Adsc and include only V because Sp and qp 
are defined on a bond, and Asdw does not include V 
because < >= < "-jo" > is independent of m. 

The pure 7r-triplet pairing state with A^^ ^ and 
Adsc — Asdw = is always energetically unfavorable 
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compared with the pure dSC state with Adsc 7^ and 
^SDW = — 0. However, since the coexistence of dSC 
and SDW {Adsc 7^ and Asdw 7^ 0) generally results 
in nonzero to {and nonzero here), the self-consistency 
of mean field calculation requires the coiisideration of tt- 
triplet pair into account from the outset.crcl' The impor- 
tant fact that the coexistence of spin-singlet Cooper-pair 
and SDW always leads to nonzero spin-triplet pair am- 
plitude with finite total momentum had been recognized 
by Psaltakis et al. in a slightly different context .EJ' The 
close relationship among the order parameters of dSC, 
SDW and 7r-triplet pair is discussed in Appendix [^. 

The mean field phase diagram in the plane of U and 
\V\ is shown in Fig. 2. While the dSC state is stabi- 
lized for small U/t, the coexistent state with dSC, SDW 
and TT-triplet pair is possible for large U/t, and the phase 
boundary between these two states is shown by solid line. 
Although the pure 7r-triplet pairing state cannot be sta- 
bilized, TT-triplet pair can condensate as a result of the 
coexistence of dSC and SDW. Since there is attractive 
interaction for spin-triplet channel in the present model, 
the coexistent region of dSC and SDW is widened by the 
inclusion of 7r-triplet pair. We note that the Fermi sur- 
face remains in the SDW state near half filling. As we 
saw in I, when SDW appears first as the temperature is 
lowered, the coexistent state can be stabilized at lower 
temperature, especially at T = 0, near half filling. 



4 
3.5 

3 
2.5 

3 2 
1.5 



n=0.9,t'/t=-1/5 

dSC+AF+Ji--triplet, 
AF 




Fig. 2. The mean field phase diagram at T = for C/ > 0, V < 0, 
n = 0.9 and t'/t = —1/5. Solid line stands for the boundary 
between the coexistent state and dSC state in the mean field 
approximation. Dotted and broken lines stand for the bound- 
ary in the right side of which the phase separation occurs in 
the random phase approximation for V" = and V' = \V\/2, re- 
spectively, where V' is the next-nearest-neighbor density-density 
interaction. N stands for the normal state. 



Generally, in the presence of finite-range attractive 
density-density interaction, the system can be hampered 
by the phase separation (PS). In order to examine the 
PS transition, we calculate the charge compressibility, 
K, given by static and uniform charge susceptibility, in 
the ground state, i.e., dSC or coexistent state. This 
phase boundary is determined from k^^ = 0. Here we 
use the random phase approximation (RPA), and take 
the RPA diagram, shown in Fig. 3, into account. The 
explicit form of k in the RPA is shown in Appendix 
B|~Jhe PS transition line is shown in Fig. 2 by dot- 



ted line. It is seen that the coexistent state is severely 
suppressed by PS but survives. This PS is expected to 
be suppressed if we take the long-range Coulomb repul- 
sion into account. Here, for simplicity, we consider the 
nexi-nearest-neighbor density-density repulsion V' > 0, 



V 



i+l ]\[ 



ii 1 
Vq —IV' cos cos qy, 



(3.4b) 



where I = ±(2; + y), ±(5; — y), and incorporate V into 
the RPA calculation, i.e., we replace Vq with Vq + V^. 
We note that this replacement, which does not alter the 
mean field equations, can bring about the nontrivial ef- 
fect on K in the coexistent state, from eq. (BTOa). For 
V' = |T^|/2, the PS transition line is shifted, as shown in 
Fig. 2, from dotted line to broken line. It is seen that the 
long-range Coulomb interaction does lead to the suppres- 
sion of the PS, which is less prominent in the coexistent 
state than in the dSC state. 




Fig. 3. The RPA diagram for charge susceptibility. 



The calculatioa of k in the RPA had been carried over 
byJVIicnas et al.^ only in the normal state. Dagotto et 
ain have shown based on quantum Monte Carlo (QMC) 
simulation that (1) PS drastically reduces the size of the 
mean field dSC region and (2) the enhancement oidx2^y2 
pairing correlation itself is not found. This QMC result 
of (2) is different from the present mean field calculation. 
We note that this QMC calculation is limited to the case 
of half filling n = 1 and relatively high temperature T = 
t/6. 

The coexistence of dSC and SDW near half filling has 
been found also in the t-J model in the slave-boson mean 
field approximationlliP and_hjL.iise of variational Monte 
Carlo calculation (VMC),E1B0 and iiiJha repulsive 
Hubbard model {V 0) by use of VMCHti In these 
studies, however, 7r-triplet pair has not been taken into 
account. The effect of 7r-triplet pair on the coexistence of 
dSC and SDW has been recently examined by Arrachea 
et alE' based on a generalized Hubbard model. In the 
repulsive Hubbard model, the nearest-neighbor hopping 
term is modified as the correlated one. 



Hrh — — 



E 

<ij>a 



{4 



tBEnian-ja 



+tAB[nia{l - Uja^) + (1 - niff)nj-^] } , (3.5) 

where three hopping integrals, tAA, tsB and tAB, incor- 
porate many-body effects into one-particle hopping pro- 
cesses phenomenologically. tAA and tsB do not change 
the number of doubly occupied sites, and tAB does, as 
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shown in Fig. 4. It is to be noted that Hch can be rewrit- 
ten as 



Hch= + c|^cj^| 



<ij>(J 



X {-t + t2{iiia + rija) + tsrii^rij-gr} , (3.6a) 



where 



t = tAA, h = tAA^tAB, h — '2tAB —tAA — tBB- (3.6b) 

The t2 terip_jC|aii_he also deduced from the bare Coulomb 
interactior£3'U'H or by including the effects of phonon 
in the antiadiabatic approximation M — > (where M 
is the phonon mass))!^) and the term describes the 
three-body interaction. Arrachea et al. have shown in 
the mean field approximation for Iab > tAA = ^bb (i-e., 
i2 < and = —2t2 > 0) and t' — that the coexistence 
of dSC and SDW is possible but prevented by 7r-triplet 
pair (and ruled out for large J7), due to repulsive spin- 
triplet pairing interaction. Hence, the effect of 7r-triplet 
pair on the coexistence of dSC and SDW is different from 
that in the present extended Hubbard model. 



i j 
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Fig. 4. The correlated hopping processes. 



3.2 Nearest- Neighbor Repulsion V > 

Next we treat the case V > 0. In this case, not only 
charge- and spin-density- wave states (CDW and SDW), 
but also orbital antiferromagnetic (OAF) and spin ne- 
matic (SN) states, in which the local staggered-cuixents 
of charge-ajid spin circulate, respectively,E3'Lil'E3'E5' are 
expected.Ej' Here we include ferromagnetism (FM) for 
the reason as we shall describe later. The order param- 
eters are 



< n.ia- >= - + erf + {p + am) cos(Qi?i), 



< c, 



-■i+p 



.a >= (.9p + <^lp) cos{QRi 



(3.7a) 



(3.7b) 



where /, p and m are real, g-p = — and l-p — —l*^. /, 
p and m describes FM, CDW and SDW, respectively. We 
take gp and Ip to be ofi-cj2.2_j^2wave symmetry, which is 
favored near half fillingES' In this case, gp and Ip become 
pure imaginary, = -g±y = ig and l±i; = -l±y = il, 
where g and I are real. The states, pith, n-^ and / ^ 
are called as OAF and SN onesjHJ'B'tSE^ respectively, 
in which there exist the staggered local currents of charge 



and spin on a bond (i, i + p)^ 

< jh±x >=- < ji.^±y > « gcos{QR,), (3.8a) 

< - < if,,±^ > cx IcosiQR,), (3.8b) 

where 



ji,i+p — ^ ^'^i'^ia'^i+P;'^ ^i+p,(7'^i'^) ' 



(3.9) 



and Vc = 1 and Vg = <J. In the OAF (SN) state, the 
local staggered current of charge (spin) circulates around 
the plaquettes, as schematically shown in Fig. 5, and 
the bond-ordered wave (BOW) does not exist, i.e., < 




Fig. 5. The OAF (SN) state in which the local staggered charge 
(spin) current circulates around the plaquettes. 



The order parameters in the mean field Hamiltonian 
are given as follows in terms of /, p, to, g and /, 

Afa/ = C//, (3.10a) 

^CDW = (8y - C/)p, ^sDW = f/TO, (3.10b) 

^OAF = 2Vg, AsN = 2Vl, (3.10c) 



where Aqaf and A 



SN 

Ip are defined on a bond 



include only V because gp and 
and ApM does not include 
V because < rii >— J2a "^^o" independent of /. 
The order parameters of CDW, SDW, OAF and SN are 
closely related to each other, which is shown in Appendix 

The mean field phase diagram in the plane of U and V 
is shown in Fig. 6. We have found that a coexistent solu- 
tion with nonzero Acdw, Asdw ApM can be sta- 
bilized for n 1. This state is ferrimagnetic, as shown in 
Fig. 7. We note that the coexistence of CDW and SDW 
{AcDW 7^ and Asdw 0) generally results in nonzero 
/ [and nonzeroAi^M here), which had been indicated by 
DzyaloshinskiicB' based on a qualitative symmetry anal- 
ysis. This is the reason why we take FM into account 
from the outset. With the present choice of parameters, 
pure FM state with only ApM cannot be stabilized, 
but FM can arise as a result of the coexistence of CDW 
and SDW. In the present case, the coexistent region of 
CDW and SDW is widened by the inclusion of FM. We 
note that the Fermi surface remains in the CDW or SDW 
state near half filling. 

It is to be noted that neither OAF nor SN can been 
stabilized solely in the mean field approximation, inde- 
pendent oft', U, V, T and n. This conclusion is contrary 
to that of Chattopadhyay et al!E3 that pure OAF or SN 
state has lower ground-state energy than pure CDW or 
SDW state for the half-filled case by introducing finite 
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t' . Moreover, a state where local-current (OAF or SN) 
and density-wave (CDW or SDW) coexist cannot be also 
stabilized. 

n=0.9, t'/t=-1/5 



U/t=4.0,V/t=1 .5,n=0.9,t'/t=-1/5 



- SDW 




Fig. 6. The mean field phase diagram at T - 
n = 0.9 and t' /t = -1/5. 



for [/ > 0, y > 0, 




(7i/2,7t/2) (0,0) 



Fig. 9. The energy dispersion relative to the Fermi level in the 
coexistent state with CDW, SDW and FM for U/t = 4.0, V/t = 
1.5, n = 0.9 and t' /t = —1/5. Full (dotted) lines stand for that 
of electrons with up (down) spin, respectively. The lower band 
of electrons with up spin is fully occupied. 



Fig. 7. The ferrimagnetic coexistent state. Each lattice site is 
shaded according to electron density. The length of each arrow 
is proportional to the magnitude of local spin moment. Lattice 
sites with larger local spin moment have higher electron density. 



For U/t = 4.0, the V dependences of Acdw, ^sdw^ 
A pM and the difference between the energy of the pure 
state (CDW for U < AV or SDW for U > W), Ep, and 
that of the coexistent state (CDW-hSDW-^FM), E^, are 
shown in Fig. 8. In the coexistent state with CDW, SDW 
and FM, IAsdh'I is larger than \Acdw\ for U > W 
and vice versa for U < 4V, and the first-order phase 
transition occurs a.t U — AV. For fixed U, |Ai^-M| rapidly 
saturates as a function oi V. It is seen that the energy 
gain in the coexistent state is very small. 

The energy dispersion in the ferrimagnetic coexistent 
state is shown in Fig. 9. There are four energy bands, 
and the Fermi surface remains as in the CDW or SDW 
state. However, the lower band of electrons with major- 
ity spin (up spin for ApM > 0) is fully occupied and the 
Fermi level crosses only the lower band of electrons with 
minority spin (down spin for ApM > 0). Therefore, this 
coexistent state is half metallic. 

The coexistence of CDW and SDW with same wave 
vectors has also been found in a ID modified Hubbard 
model ftp^a quarter-filled band in the mean field approx- 
imation£a) In this coexistent state, the wave vector of 
charge and spin density, and the phase difference be- 
tween CDW and SDW, AO, are equal to 2kF = n/2 and 
7r/2, respectively, and the magnitude of local spin mo- 
ment is equal at each site, as shown in Fig. 10. On the 
other hand, in our coexistent state, q = Q = (tt, tt) and 



A6 = 0, and the magnitude of local spin moment is dif- 
ferent at each sublattice, as shown in Fig. 7. This ferri- 
magnetic coexistent state is the 2D version of that found 
in the 3D Hubbard model {V — 0) which hacLbeen de- 
noted as the special ferrimagnetic (S.F.) state.EZP In the 
present 2D case, this coexistent state ioi V = can be 
stabilized for 12<,C//t<,14 (not shown in Fig. 6). 



Fig. 10. The coexistent state with 2kp CDW and 2tH SDW 
found in a quarter-filled ID modified Hubbard model.Cj' Each 
lattice site is shaded according to electron density. The length of 
each arrow, proportional to the magnitude of local spin moment, 
is equal at each site. 



§4. Renormalization Group Analysis 

In the last section, we examined possible ordered states 
for U > and V ^ in the mean field approxima- 
tion. In this section, we examine the effects of fluctu- 
ation on these ordered states which are not taken into 
account in the mean field calculation. As a theoret- 
ical treatment beyond the mean field level, we adopt 
the renormalization group (RG) method for the saddle 
points whirl-i "has been applied to the Hubbard model 
{V = O),^^^^^^'© and determine the most dom- 
inant correlation in the normal state. We also discuss the 
possibility of the coexistent states beyond the mean field 
approximation. 

4-1 Saddle Point Singularity 

We consider the special case where the Fermi level in 
the absence of interaction lies just on the saddle points 
Qa = (7r,0) and Qs = (0,7r), i.e., /i = 4t' (n - 0.83), 
and focus on electrons at these two saddle points on the 
Fermi surface, just as two Fermi points in ID electron 
systems. The Fermi surface is shown in Fig. 11. 

First, we examine the behavior of the following 
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U/t=4.0, n=0.9, t'/t=-1/5 U/t=4.0, n=0.9, t'/t=-1/5 




V/t V/t 
(a) (b) 



Fig. 8. The V dependences of (a) Aq£)\y, '^sdw ^fm ^^rid (b)the energy difference Ep — Ec (scaled by t) between the pure state 
(CDW or SDW) and the coexistent state (CDW+SDW+FM) for U/t = 4.0, n = 0.9 and t' /t = -1/5. 



t'/t=-1/5, Mf 




logarithmically divergent 



C , 2 -Sc „ C 



K9 



-f-n log — for w < rEc 



-Pi 



for LO » rEc 



-r^ log £^ for w < rEc, 
-Ki for uj > rEc-, 



(4.5a) 
(4.5b) 

(4.5c) 



Fig. 11. The Fermi surface in the absence of interaction for t' /t ■ 
— 1/5 and = 4t', in which case n ~ 0.83. 



particle-particle (K) and particle-hole (P) correlation 
functions, 

ifac' = hm / Gc{k,e)Ga'{-k + q,-e), (4.1a) 



where i?c > and > (w <C i?c) are the ultraviolet 
and infrared energy cutoff, respectively, 

(4.6a) 
(4.6b) 



i + \/r 



log 
1 

2^ 



2r 

arctan( 



2r 



Vl - 4r2 ' 



(4.6c) 



lim 



for a, a' = A, i3, where 
G„(fc,e) 



/ Ga{k, e)Ga'{k + q, e), (4.1b) and r = \t'\/t. c, c' and c" as a function of r are shown 
^fc,e in Fip, 12_ F,snpr,iallv for small t-_ 



in Fig. 12. Especially for small r, 

c, c" ~ 1, c' ^ — logr. 



(4.7) 



1 



is the one-particle Green function in the absence of in- 
teraction for electrons near the saddle point Qa and 

L^Xk.. (2^/2^' ^^'^^ 

where the cutoff around the saddle points, kc, is intro- 
duced. We note that 



K,=Kaa = Kbb, K2 = K 



AB 



K 



BA, 



(4.4a) 



stand for Cooper- and 7r-pair correlation, respectively, 
and 



Pi = Paa = Pi 



BB, 



P9,^P^ 



AB 



P, 



BA, 



(4.4b) 



for uniform and staggered density-density correlation, re- 
spectively. For /i — At', these correlation functions are 



For r > Tc 0.276, c > c' and Pi is more divergent 
than P2 for a; — > 0. We note c" <max{c, c'}, i.e., 7r-pair 
susceptibility K2 is always less divergent than particle- 
hole susceptibility. For t' /t = —0.2, c" < c < c' and they 
are comparable in magnitude. 



4-2 Renormalization Group Method for the Saddle 
Points 

In the last subsection, we saw that the saddle points 
on the Fermi surface lead to logarithmic divergence of 
particle-particle and particle-hole correlation functions. 
This implies that the fluctuation effect becomes strong. 
In the RG approach, we assume that the single renormal- 
ization group variable x = log determine the behavior 
of the system. The increase of x represents renormal- 
ization towards lower energy scale. For simplicity, we 
neglect (1) the deformation of the Fermi surface by in- 
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Fig. 12. c, c' and c" as a function of r. 



teraction and (2) fc-dependence of interaction for small 
\k\ < kc, i.e., we consider only eight coupling constants, 
gis {i — 1,2,3,4 and s =_L, ||). This interaction of the 
g-61ogy type is shown in Fig. 13. gi and g^ (52 and g^) 
stand for the backward (forward) scattering processes 
with large (small) momentum transfer, respectively. Es- 
pecially, gi and gs describe 'exchange' and 'Umklapp' 
processes. In I, only gi_i and g3± were treated and taken 
to be momentum-independent all over the magnetic Bril- 
louin zone™ 



gl 



S3 



S2 



g4 



Fig. 13. The scattering processes. SoUd and dashed lines stand 
for electrons near Qj^ = {tt,0) and Qg = (0, tt), respectively. 



We start with the renormalization of the couplings 
in the one-loop approximation. One-loop diagrams are 
shown in Fig. 14. The scaling equations are 



911- 
9i\\ 

92\\ 

5311 



-2.gi_L(.gi|| - .g2||)A, 
'2.9111,9211^2 - 2.9i||54||Pi 

-[25111(5111 - 92\\) + {gf± - gl\\ 

+ (53\-53||)]A, 



(4.8a) 



(4.8b) 



-(5?± + 5D-f^2 - 254_l(5i|| - 52||)A 
-{9i±+9l±)P2, (4.8c) 
-(5?|| +52||)^2 - 2(.g4||.gi[| - 54±.92_l)A 
-(5211+ 5311 ) A, (4.8d) 
-253_L54_L-^'i - 253_l(52_l + 52[| - 5il|)^2(,4-8e) 
-2g3||54||ifi - 2(253||52|| - .93±5i±)A, (4.8f) 



54-L 



5411 



igL 



'gi±)Ki 



^[gli 



gli 



252±(5i|| -52|l)]A, 



(4.8g) 



-54||)^i 
(25411 



94A 



25211(5111 -5211) 



-(52\-52||)]^i' 



(4.8h) 



which are to be solved with the initial conditions gis {x — 
Xi) = g'is {' =d/da;), where are the bare coupling 
constants. For example, g^x^ has the following form to 
one loop order, 

g3± = -2.93^34% A'l - 233^(52% + - 5?|| )^2. (4.9) 

By differentiating this equation by x and replace g^, by 
5is, i.e., the bare coupling constants by the renormalized 
ones, we obtain the scaling equation eq. (4.8e). 
If we take g^^ = as the initial conditions, the relation 
gi± = gi\\ holds all through the flow. Therefore, the 
above scaling equations are simplified as follows,^ 



51 

52 

53 

54 



-2.9i.92i^2 - 2gigiPi - 2.91(32 - .9i)P2(4.10a) 
-(5? +51)^2 + 254(52 -5i)A 
'{9I ^ • 



.93)^2, 
-2.93.94X1 - 2.93(2^2 
-{gl+gl)K^ 
-[.9? + .94 + -2.92(52 



5i)A 



51)] A. 



(4.10b) 
(4.10c) 

(4.10d) 



where = gi±_ = 5i||- 

The divergence of gis (x) at a finite x indicates the exis- 
tence of the strong coupling fixed point, i.e., signals the 
development of an ordered state, at finite energy scale 
or finite temperature. (Strictly speaking, this finite on- 
set temperature is an artifact of the present approxima- 
tion in the 2D systems, and should be interpreted as a 
crossover temperature, or a critical temperature when fi- 
nite three-dimensionality is assumed.) The properties of 
this strong coupling fixed point can be obtained quali- 
tatively from various response functions. The response 
functions in the one-loop approximation are obtained 
from one- loop diagrams shown in Fig. 15. The response 
function, 



Ri, 



dre 



IOt 



1 



(4.11) 



/o N 

where v stands for the kind of correlation [v =dSC, 
SDW, • • •), has the following form to one-loop order, 

1 



:gARv) 



(4.12) 



where g'^ is the coupling constant (linear combination of 
(?£,) and i?° is the simple bubble. If we differentiate this 
equation by x and replace g^, R^ by the renormalized 
ones, gi, and R^, we obtain 



R,=RU1 + -g,R, 



(4.13) 



This equation is to be solved with the initial condition 
that R^{x = Xi) ^ 0. Since Ri, is positive, i?^ can 
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be divergent for > and are suppressed to zero for 
gi^ < 0. In this paper, we consider the response functions 
shown in Fig 16, where Ossc, Or, and Ops stand for s- 
wave Cooper-pair, ry-singlet Jiair with a total momentum 
Q and total spin S = 0J2j'l3 and uniform charge den- 
sity, respectively. The most divergent Rps is interpreted 
to describe the phase separation (PS). The relationship 
among these order parameters is discussed in Appendix 
[~A| . We note that each correlation is treated indepen- 
dently in the above procedure. Therefore, we can de- 
termine the most dominant susceptibility in the normal 
state, and cannot assess the coexistence of different or- 
ders. 




Fig. 14. Diagrams contributing to the one-loop order correction 
to coupling constants. 




(a) 




Fig. 15. Diagrams contributing to the one-loop order correction 
to response functions. 



4-. 3 Phase Diagram 

We solve the scaling equations, eq. ( |4.8| ) and (4.13) 
with the initial conditions, 



5?^ = 



.93% = U + 2Vq = U~AV, 
.g°j^ = U + 2Vo^U + W, 
. gO,l = 2Vq = -4y, 



2Vq = 41/, 



(4.14a) 
(4.14b) 
(4.14c) 
(4.14d) 



at a; = a:^. Here, we take Xi = Q and Ri,{xi) = for 
simplicity, although the solution of the scaling equations 
depends on the value of Xi and R^{xi). 

Before we show our results, we refer to previous results 
ioT U > and V = obtained by many authors. Hjj 



V 








UOV_y 


1 




93L — g4± 


sSC 


1 


4Ki 


-(53^ +54^) 


7r 




4Jf2 


Sll - S2± 






4:1(2 


-(Sl± + S2±) 


CDW 


Ep<r (y^Cp+Qa 


-4P2 


-(Sii + S3±) - S 


SDW 


Spa- ^Cj,^Cp_|_go- 


-4P2 


Si± -1- S3± - s 


OAF 


i Spa ffpcJo-Cp+Qo" 


-4P2 


-(Sii - S3±) - s 


SN 




-4P2 


(9i± - S3±) - S 


FM 


X^PCT ^CpffCpCT 


-4Pi 


S2± + S4± + S 


PS 




-4Pi 


-(S21 -1- S4±) -1- S 



Fig. 16. Operator Ov, bare response function R'l{> 0) and cou- 
pling constant gi, for each correlation u. Each response function 
can be divergent (or there exists mean-field solution < Oi, 0) 
for gu > 0. g is defined as g = giy — g2\\. In the RG method 
where only electrons near the saddle points are taken into ac- 
count, the sum over p is restricted to p ~ (""lO), (0, vr) and the 
p-dependence of the d^2_„2-wave factor Wp oc cospi, — cos py is 



ignored. In our calculation, we take Wp = sgn{cospi, 



iPy). 



can be rewritten as follows, 

— n^m 



H, 



U 



u 



2 ^' 



(4.15) 

If we regard the second term in the r.h.s of eq. ( 4.15| ) 
the chemical potential shift, we can take g'^ 



9^\\ 



as 
U 



as the initial conditions and therefore use eq. ( 4.10| ) 
as the scaling equations of the coupling constants. For 
the perfect nesting case r = \t'\/t — 0, SchulzEEP and 
DzyaloshinskiiEB' showed that SDW occurs, and pointed 
out that smalL deviations from halLiilling lead to dSC. 
Lederer et alS3 and Furukawa et alS3> solved the scaling 
equations foc-t <c; 1 and the same results. Especially, Fu- 
rukawa et alS3 indicated that the correlation of 7r-triplet 
pair is suppressed to zero. In these calculations, however. 



Pi and K2 are neglected in eq. (I.IC) for the reason that 
they are less singular than P2, i.e., c, c" ^ c' for r <C 1, 
iupea. ( |4.5| ) and ( |4.6| ). On the other hand, Alvarez et 
al£3 solved the flow equations with the initial condition 
g^j_ = U and = 0, by neglecting the generation of 
gfjii and omitting particle-particle diagram Ki and K2 in 
eq. (Q, and showed that dSC, SDW and FM can be 
stabilized depending on U/t and r. 

Now, we show the phase diagram in the plane ofU > 
and V for t' = —1/5 in Fig. 17, in which the ordered 
state with highest onset temperature is shown. First, 
we discuss the case that each ordered state is treated 
independently, because we cannot examine the coexis- 
tence of different orders in the RG method. The crucial 
difference from our mean field result is that supercon- 
ductivity appears even for [/ > and V > 0, i.e., dSC 
for U > AV and sSC for U < AV. This result that dSC 
is possible for small > as well as V = near half 
filling is consistent with a recent calculation based on 
the fluctuation-exchange (FLEX) approximation.E3' Ex- 
cept for superconductivity for U > and > 0, the RG 
phase diagram is qualitatively same as the mean field one 
when we do not take the coexistence of different orders 
into account, i.e., SDW and CDW appear for U > AV 
and I] < AV, respectively, and attractive V < favors 
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dSC for small \ V\ and PS for large \ V\, respectively. With 
regard to the correlation of 7r-triplet pair, our RG calcu- 
lation has shown that it can be divergent for attractive 
y < and large \V\ but is always subdominant. Sim- 
ilarly, FM cannot be the most dominant solely. These 
results are also consistent with our mean field ones. 

Next, we discuss the possibility of the coexistence of 
different orders at low temperature, especially at T = 0. 
It is very important that our RG calculation shows the 
existence of a region where the onset temperature of 
SDW or CDW becomes highest, as shown i n F ig. 17. 
Since our mean field calculation in §3.1 or §3.2 shows 



that the Fermi surface remains in the SDW and CDW 
states near half filling, we might expect to find a sec- 
ond phase transition at lower temperature in such SDW 
and CDW states. Therefore, at lower temperature in 
the SDW region in Fig. 17, (1) the coexistent state with 
dSC, SDW and 7r-triplet pair found for V < in the 
mean field approximation might be expected to survive 
for not only V < but also > 0, and (2) the ferrimag- 
netic coexistent state with CDW, SDW and FM found 
for U > AV > in the mean field approximation might 
be expected to survive. In fact, as we have pointed out 
in §3.1, VMC calculations for [/ > and V = shojssi.th£ 
coexistence of dSC and SDW at low temperature Jl3'L§ 
although TT-triplet pair has been neglected. Similarly, at 
lower temperature in the CDW region in Fig. 17, (3) the 
ferrimagnetic coexistent state with CDW, SDW and FM 
found for AV > C/ > in the mean field approximation 
might be expected to survive. Especially, 7r-triplet pair 
(FM), which cannot be stabilized solely in the parame- 
ter region considered in the present mean field and RG 
approximation, might be expected to arise as a result of 
the coexistence of dSC and SDW (CDW and SDW). In 
order to assess the effect of 7r-triplet pair (FM) on the 
coexistence of dSC and SDW (CDW and SDW), we need 
another theoretical treatment. 



4 

3.5 
3 

2.5 
2 

1.5 
1 

0.5 




tVt=-1/5 









SDW / 


\ dSC 




CDW 


PS \ 






/( A 






/ sSC 





-1.5 



-0.5 




V/t 



0.5 1 



1.5 



Fig. 17. The RG phase diagram for t' /t = -1/5 and ^ = At' 



Here, we refer to the possibility of the coexistence of 
sSC and CDW. In the g-61ogy model in I, it can be eas- 
ily shown that the coexistent state with sSC, CDW and 
?7-singlet pair can be stabilized nearJialf filling at low 
temperature for < and gn_ < O.Ef As seen from eq. 



in the present extended Hubbard model. Therefore, this 
coexistent state might be expected to be stabilized also 
in the extended Hubbard model for [/ < and y > in 
the mean field approximation. Moreover, based on the 
above discussion, it might be expected to survive in the 
presence of fluctuation for not only U < but also U > 0, 
at lower temperature in the CDW region in Fig. 17. This 
will be reported elsewhere. 

Finally, we refer to the ambiguities of the above RG 
method. Since there exist not only log- but also log^- 
divergence in the particl e-p article and particle-hole cor- 
relation functions, eq. (4.5), we cannot safely take the 
limit — > in the scaling equati ons of cou pling con- 

and (1.13|), i.e.. 



slants and response functions, eq. (4 



it is not clear at all whether the above RG treatment is 



(2.2), this case corresponds to that of U < and V > 



valid or not. In fact, the solution of eq. (4^) and ( 4. IS ) 
depends on the initial value Xi . If we consider only the 
most singular log^ term in Ki (and P2 for r = 01 and 
take y = — log^ as a new scaling variablej23^ we 
can safely take — > limit in the scaling equations of 
coupling constants. In this case, the above RG method 
might correspond to a parquet summation of leading J 
divergences, rather than renormalization procedure.! 

§5. Conclusion and Discussion 

We have studied in detail possible ordered states, es- 
pecially possible coexistence of different orders, near half 
filling in the 2D extended Hubbard model with on-site re- 
pulsion U > and nearest-neighbor interaction V, with 
emphasis on electrons around the saddle points (7r,0) and 
(0,^). 

First, we have determined the phase diagram at T = 
in the mean field approximation. For V < 0, we have 
shown that the coexistent state with dSC, SDW and n- 
triplet pair can be stabilized near half filling. Here, we 
have indicated the following important fact which has 
often been neglected in previous studies: when we dis- 
cuss the coexistence of dSC and SDW, it is necessary 
to take TT-triplet pair into account from the outset, be- 
cause in general the coexistence of dSC and SDW results 
in TT-triplet pair and is affected by TT-triplet pairxi Espe- 
cially, TT-triplet pair, which cannot condensate solely in 
the present model, can arise through the coexistence of 
dSC and SDW. Since the phase separation (PS) is gen- 
erally expected to occur in the presence of finite-range 
attractive interaction such as y < 0, we have examined 
the effect of PS on the mean field ground state in the ran- 
dom phase approximation (RPA). The coexistent state 
with dSC, SDW and 7r-triplet pair is severely hampered 
by PS but survives, and that the long-range Coulomb 
repulsion such as next-nearest-neighbor density-density 
repulsion leads to the suppression of PS. On the other 
hand, for V > 0, we showed that a ferrimagnetic co- 
existent state with commensurate charge-density-wave 
(CDW), SDW and ferromagnetism (FM) can be stabi- 
lized near half filling. Here, we have indicated the fol- 
lowing important fact: when we discuss the coexistence 
of CDW and SDW, it is necessary to take FM into ac- 
count from the outset, because in general the coexistence 
of CDW and SDW results in FM and is affected by FM. 
Especially, FM, which cannot be stabilized solely with 
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the present choice of parameters, can arise through the 
coexistence of CDW and SDW. It is to be noted that the 
above mean field coexistent states near half filling can be 
stabilized at low temperature, especially at T = 0, when 
CDW or SDW, in which the Fermi surface remains near 
half filling, arises first at high temperature. 

In order to examine the effects of fluctuation on the 
mean field ordered states, we have adopted the RG 
method for the special case that the Fermi level lies just 
on the saddle points. We have shown that the crucial 
difference from our mean field result is that supercon- 
ductivity can arise even for U > and V > 0; dSC 
and sSC for U > W and U < iV, respectively. Except 
for this difference, the RG phase diagram is qualitatively 
same as the mean field one when we do not take the co- 
existence of different orders into account, e.g., SDW and 
CDW can arise for U > AV and U < 4V, respectively. 
Especially, the correlation of 7r-triplet pair or EM cannot 
be the most dominant solely. Here, it is very important 
that a region where the onset temperature of SDW or 
CDW becomes highest is found in the RG phase dia- 
gram. Since the Fermi surface remains near half filling 
in these SDW and CDW states, we might expect to find a 
second phase transition at lower temperature. In the RG 
method, however, we cannot assess such possibilities. On 
the other hand, the mean field approximation, which is 
often questionable for the 2D case, is of great advantage 
in that we can study the stability of coexistent states 
with different order parameters quantitatively. There- 
fore, we can conclude that our mean field calculations 
indicate the possibilities that (1) SDW, which has been 
shown in our RG calculation to arise first at high temper- 
ature for U > 4:V, coexists with dSC and vr-triplct pair, 
or with CDW and FM, at lower temperature, and that 
(2) CDW, which has been shown in our RG calculation 
to arise first at high temperature for U < AV, coexists 
with SDW and FM, at lower temperature. 

Throughout this letter, we have assumed YBCO-type 
Fermi surface by introducing t' and consider only com- 
mensurate (C) SDW or CDW. Near half filling, however, 
incommensurate (IC) ordering or stripe formation can be 
expectedrJispecially for t' = in the repulsive Hubbard 
modelB3'E3E3) Recently, the effect of V on such stripe 
states has been examined. Eil* The effect of IC ordering 
on the stability of the coexistent states (with dSC, C- 
SDW and 7r-triplet pair, and with C-CDW, C-SDW and 
EM) is beyond the scope of the present study. In the xe- 
pulsive Hubbard model with t' = 0, Giamarchi et alc3 
have shown that the coexistent state with dSC and C- 
SDW state have higher energy than that with dSC and 
IC-SDW. 

With regard to dSC, we have assumed that the su- 
perconducitng gap symmetry is purely of da;2_y2-wave. 
However, there are a few indications that dSC mixed 
with pj^pnmonpnts of other symmetry can be stabi- 
lized,Bc3'Oc3 dependent on interaction, electron den- 
sity, etc. Such mixed pairing states leave much room for 
future studies. 

The author thanks to H. Fukuyama, H. Kohno and M. 
Ogata for valuable discussions. 



Appendix A: Symmetry among Various Order 
Parameters 

We examine the relationship among the order param- 
eters as shown in Fig. 16. First, the operators of dSC, 
SDW and 7r-triplet pair are equivalent in that they are 
'rotated' into each other. 



[6dsc,dl]^2dsDW- 



(A-lb) 



This relationship underlies the SO{5) theory, in which 
dSC and SDW are unified into a /iwe-dimensional vec- 
tor superspin and describes excitation towards the 
SDW (dSC) directio n in th e dSC (SDW) ground statefi) 
We note that eq. ( A-lb| ) holds due to — 1, i.e., 

Wp = Sgn(c0SPa; — COS Py). If Wp OC COSPx — COS Py, it 

is satisfied only in the long wavelength limit where only 
electrons near the saddle points are important. It is very 
important to note that if two of the above three order 
parameters coexist, another one also results generally. 
As we have already shown in the mean field approxima- 
tion, the coexistent state with dSC, SDW and 7r-triplet 
pair can be stabilized near half filling in the 2D extended 
Hubbard model for C/ > and V < 0. 

Next, the operators of sSC, CDW and 77-singlet pair 
are equivalent in that they are 'rotated' into each other, 



[OcDW, disc] =261 [6^,OcDw] = 2d,sc, (A-2a) 



sSC- 



20, 



CDW, 



(A-2b) 



This relationship underlies the SO{3) theory, in which 
sSC and CDW are unified into a t/iree-dimensional vec- 
tor pseudospin and Orj describes excitation towards 
the CDW (sSC) direction in the sSC (CDW) ground 
state. oEj^ In fact, such properties are useful in the at- 
tractive Hubbard model. It is very important to note 
that if two of the above three order parameters coexist, 
another one also results generally. It can be easily shown 
in the mean field approximation that a state where sSC, 
CDW and 77-singlet pair coexist can be stabilized near 
half filling iii-the g-61ogy model used in I for 173^ < 
and gi± < OB 

We note that similar relationships hold among dSC, 
OAF and 77-singlet pair, 

[6lddsc]=2idoAF, [doAF,dr,]^2iddsc, (A-3a) 



rfSC OoAF 



2iOl 



(A-3b) 



and among sSC, SN and 7r-triplet pair. 



[dlsc,0.] = 2idsN, [d,N,d,sc] ^2id^, (A-4a) 



SN\ 



2iO 



t 

sSC- 



(A-4b) 



If Wp cx cospx — cospj,, eq. (A-3b) and (A-4b) hold 
only approximately. In each case, if two of the above 
three order parameters coexist, another one also results 
generally. It can be easily shown in the mean field ap- 
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proximation that a state where dSC, OAF and 77-smglet 
pah (sSC, SN and 7r-triplet pair) coexist can be stabi- 
hzed in the g-61ogy model used in I for g^^ > and 
9i± < (53i < and gii_ > 0)B 

Moreover, eq. ( |A-lD -( A-4) imply close relationships (1) 
between sSC and dSC, (2) between ry-singlet pair and tt- 
triplet pair, and (3) between density-wave (CDW, SDW) 
and local-current (OAF, SN). In fact, the hermite opera- 
tors of uniform bond-charge and bond-spin (BC and BS) 
density with (i2;2_j,2-symmetry defined as 



Obc = X! '^p^U^pcr, Obs 

pa 



satisfy similar commutation relations, 



(A-5) 



(A-6a) 



[dfi,d^]^2da, (A-6b) 
for (a,/3,7) =(sSC, dSC, BC) and (77, tt, BS), and 

[d^,dp\=2id^, [6^,6^] =2idp, (A-7a) 

[6^,0^] = 2i6„ (A- 7b) 

for (a,/3,7) =(SDW, BC, SN), (CDW, BC, OAF), 
(SDW, BS, OAF) and ( CDW BS, SN). If Wp cx cosp^^ - 
cospj,, eq. ( A-6b ) and ( A-7b| ) hold only approximately 
in each case. 

The relationship among order parameters which are 
rotated into each other by O,,, Ojr, Obc and Obs is 
summarized as shown in Fig. 18, where the connection 
with Odsc and Ossc by Obc is not explicitly shown. 

dSC 

71 / \ ri 



SDW< 

A 

BC 

SN < 
y 



BS 



BS 



> OAF 

A 

BC 

► CDW 

-i 



71 T] 

sSC 



Fig. 18. The relationship among various order parameters. 



Appendix B: Charge Susceptibility in the Ran- 
dom Phase Approximation 

The charge compressibility k is equal to static and 
uniform charge susceptibility Xi 



K = lim X iq,q,UJ 



0), 



(B-1) 



where x^ili q' t^) is the retarded density-density corre- 
lation function, which is obtained from the analytic con- 



tinuation of xili q' 1 i"-^;) through iuji lo + 16, 

. 

X{q,q\icji)^ dTe'"'^x(9,?',^)- (B-2a) 
Jo 

x{Qi1'tT) is the two-particle thermal Green function in 
imaginary time given by 



dq, q', 't) = X! 



-q'a' 



(B-3a) 



>, (B-3b) 



where Uqa = Y.kAa'^k+qa, n^g^ = {Uqa, 

Ncra' in the matrix form as follows. 



N 



^ We write 



(B-4) 



For the mean field Hamiltonian with Adsc, Asdvk and 
AjT, the one-particle thermal Green functions, defined by 

G^^'(p,p',t) = - < TrCp^{T)cl,^, >, (B-5a) 
Fl,{P,p\r) ^ - < Trcl{T)cl,^, >, (B.5b) 

Faa'iPiP' ^'t) = - < TrCpa{T)Cp: c,' >, (B-5c) 

have the following form, 

Gcrcr{p,P,T) = 5p'^pGaal{p,T) + 5pi ^p+QGaa2{p, t) , (B-6a) 

fI^{p,p', t) ^ Sp,,^pFl^,{p, r) + Sp,.^p+QFig{p, r). (B-Gb) 
Therefore, it is seen that x has the following form, 

N{q,q' ,t) = Sq',qNi{q,r) + Sq,^q+QN2iq,T), (B-7a) 



No - 



2|| 
A^2i 



^2 



2|| 



(B-7b) 



(B-7c) 



The subscript 1 and 2 represent normal and Umklapp 
part, respectively. It is to be noted that there exists 
Umklapp part N2 for Asdw 7^ or A^r ^ 0. With 
regard to x(9, t), it has only normal part, 

x{q,q',T) = Sq'^q ■ xo{q,T), (B-8a) 

Xo{q,T) = 2{N^iq,T)+NMq,r)} • (B-8b) 

Next, we treat the effect of interaction in the random 
phase approximation (RPA). The RPA equation is dia- 
grammatically shown in Fig. 19. Here, we consider only 
(a) for interaction vertex in Fig. 19, which leads to the 
diagram shown in Fig. 3. In the frequency space, the 
RPA equation is written as follows {z = iioi), 

NRPA{q,q',z) = N{q,q',z) 

+ J2 ^(9. qi,z)9{qi)NRPAiqi,q', z), (B-9a) 
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(B-9b) 
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g\\{q) = -2V,, g^{q) = -{U + 2Vq), 



(B-9c) 



and Vq — V{cosqx + cos^j^). Since it can be seen that 
Nrpa has the same form as eq. (B-7), xrpa is obtained 
as follows, 



gtYq' 



XRPA{q,q + Q) = 0, 



(B-lOa) 
(B-lOb) 



where 



9\\iq)±9±iq), 



^ 9q+QUq U^+Q, 



(B-lOd) 
(B-lOc) 



and the variable z is not explicitly written. Therefore, 
we can obtain the charge compressibility n in the RPA, 
from eq. (BT) and (B-lOa). In the normal or dSC state, 



Xrpa has a familiar form due to = 0, 



XRPA(q,q) = -rr+ 

Xn 



Xa{q) 



l + ^{U + AVq)xoiq)' 



(B-11) 



where xo is given by eq. (B-8b). However, since ^ 
in the SDW or 7r-triplet pairing state, especially in the 
coexistent state with dSC, SDW and 7r-triplet pair, xrpa 
has a nontrivial form. 





+ 



(a) 



Fig. 19. The RPA equation for charge susceptibility. The square 
stands for interaction vertex. In our calculation, the diagram of 
(b) for interaction vertex is neglected. 
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